BLOCK-MODIFIED WISHART MATRICES AND FREE POISSON 

LAWS 

TEODOR BANICA AND ION NECHITA 

Abstract. We study the random matrices of type W — (id ® ip)W, where W is a 
complex Wishart matrix of parameters (dn,dm), and (p : M n (C) M„(C) is a self- 
adjoint linear map. We prove that, under suitable assumptions, we have the d — > oo 
eigenvalue distribution formula SmW ~ 7r m „ p M is, where p is the law of ip, viewed as a 
square matrix, it is the free Poisson law, v is the law of D = ip(l), and <5 = tr(D). 



Introduction 

A complex Wishart matrix of parameters (dn, dm) is a random matrix of type W = 
■^GG*, where G is a dn x dm matrix with independent complex A/"(0, 1) entries. In 
the limit d — >■ oo, the eigenvalue distribution of W converges to a certain law 7r t , with 
t = m/n, computed by Marchenko and Pastur in [16]. Later on, the free probability 
theory of Voiculescu [21], [22], [21] has led to a new, conceptual point of view on this 
result: the limiting law ir t is in fact the "free Poisson law" of parameter t. 

These fundamental results have been subject to a number of extensions and general- 
izations. In particular, three types of "block-modified" versions of the Marchenko-Pastur 
theorem appeared in the recent random matrix and free probability literature: 

(1) The matrices of type (1 ® E)W, where E is the diagonal matrix formed by the 
n- roots of unity, were investigated some time ago in [3]. The limiting laws here 
are certain compound free Poisson laws, called free Bessel laws. 

(2) The matrices of type {id® tr(.)l)W, where tr is the normalized trace of the nxn 
matrices, appeared in connection with the quantum information theory problems 
investigated in [TT]. The limiting laws here are the free Poisson laws. 

(3) The matrices of type (id<g>t)W, where t is the transposition, were investigated by 
Aubrun in [2]. His computation, leading to shifted semicircles, was extended in 
[3], where the limiting law was shown to be a difference of free Poisson laws. 

The main motivation for the above results comes from quantum information theory. 
In quantum information theory, the partial transposition map is known to be an "en- 
tanglement witness": it allows to test if a quantum state (represented by a positive, 
unit trace matrix) is entangled, in the following sense. If a bipartite quantum state 
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p G Mrf(C) £§> M n (C) is separable (i.e. it can be written as a convex combination of 
product states p- 1 ^ <8> pf^), then its partial transposition p r = (id ®t)p is also a quantum 
state. However, if p is entangled, then p r may fail to be positive. In the case where p r 
is a positive matrix, the quantum state p is said to be PPT (Positive Partial Transpose). 
Hence, separable states are always PPT and non-PPT states are necessarily entangled. 
The equivalence of entanglement and non-PPT is known to hold only for total dimension 
smaller than 6 (2 x 2 or 2 x 3 product systems) and it fails for larger dimensions, in 
the sense that there exist PPT entangled states. In the same spirit as in [2], the results 
in [5] regarding the positivity of the support of the limit measure can be interpreted as 
results about typicality of PPT states for large quantum systems. The Wishart matrices 
(normalized to have unit trace) are known to be physically reasonable models for random 
quantum states on a tensor product C d ®C n , the parameter m of the Wishart distribution 
being related to the size of some environment <C dm needed to define the state. So, as a 
conclusion, the various technical results in [2] and [5], not to be detailed here, indicate 
that when m > 2 and n < m/4 + 1/m, a typical state in C d ® C n is PPT. 

The starting point of the present work was the following observation: the matrices W, 
(1 eg) E)W, (id <g) tr(.)l)W and (id <8> t)W appearing in the above considerations are all 
particular cases of matrices of the form W = (id ® <p)W, where <p : M n (C) — > M n (C) 
is a linear map. So, a natural problem is to try to compute the asymptotic eigenvalue 
distribution, with d — » oo, of such general "block-modified" Wishart matrices W. 

In this paper we solve this problem, under suitable assumptions on ip. Let us first recall 
that the linear maps of type ip : M n (C) — > M n (C) are in correspondence with the matrices 
A G M n (C) £§> M n (C). We will use the following formula for this correspondence: 

(p(A) = (Tr g> id)[(t ® id) A ■ (A <g> 1)] 

Here Tr and t denote the usual trace and transposition of the n x n matrices. This 
correspondence, when restricted to the subclasses of completely positive maps, is known 
in quantum information theory as the Choi-Jamiolkowski isomorphism. See [7]. 

Given a real positive measure pi, not necessarily of mass 1, we denote by it^ the corre- 
sponding compound free Poisson law. Also, we write W ~ \x if the d — > oo asymptotic 
eigenvalue distribution of W is /i. With these notations, our main result is as follows: 

Theorem. Let W = (id <g> tp)W , where W is a complex Wishart matrix of parameters 
(dn,dm), and where (p : M n (C) — > M n (C) is a self-adjoint linear map, coming from a 
matrix A G M n (C) (g> M n (C). Then, under suitable "planar" assumptions on <p, we have 
5mW ~ Tfmnp v, with p = law(A), v = law(D), 5 = tr(D), where D = <p(l). 

This result generalizes the above-mentioned computations in [2], [5], [11], [16]. The 
proof is quite standard, first by applying the Wick formula, then by letting d — > oo and 
by using Biane's bijection in [8], and finally by using Speicher's free cumulants [T9] . 

We should mention that the combinatorics of some very general "modified Wishart" 
matrices was already investigated some time ago by Graczyk, Letac and Massam [T2] . 
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Our assumption that the modification is performed blockwise is of course a key one: this 
makes the whole combinatorics controllable, and leads to the above result. 

Let us also mention that, as kindly pointed out to us by the referee, there should be as 
well a direct proof for the above result. Indeed, let be the matrix units of the inclusion 
M n (C) C M n d(C). It is known from Voiculescu's work [23] that is asymptotically from 
W. Now the matrix W can be viewed as a linear combination of matrices of type e^Weki, 
so, in principle, combinatorial computations of free probabilistic flavour, in the lines of 
the book by Nica-Speicher [17] for example, should probably do as well the job. Several 
proofs of this type, for related results, have been worked out, e.g. in [I], [9], [15], [IB] . 

We assume in the above statement that tp is self-adjoint, in the sense that it maps self- 
adjoints to self-adjoints. It would be interesting to understand as well the non-self-adjoint 
case, as to cover the "free Bessel" computation in [3J, where <f(A) = EA, with E being 
the diagonal matrix formed by the n-roots of unity. We should mention that our formula 
5mW ~ Timnp v covers as well this situation, but only in the formal sense of [3J . 

Some other problems concern a possible relation with the complex reflection groups. 
Indeed, the exact planarity assumptions needed on <p, not to be detailed here, concern 
certain sets of partitions P even (2k,2k), for which we refer to Definition 4.4 below, and 
these sets are known to span the centralizer spaces for the hyperoctahedral group [I]. 

The above questions are probably all related. A common answer to them should prob- 
ably come from a very general random matrix/free probability formula, extending our 
present formula 5mW ~ vr mnp Kl v. But we have no further results here. 

The paper is organized as follows: in 1-2 we perform a joint study of the compound 
free Poisson laws and of the block-modified Wishart matrices, and in 3-4 we develop a 
number of supplementary ingredients, and we state and prove our main result. 

Acknowledgements. We would like to thank Mireille Capitaine, Benoit Collins, Maxime 
Fevrier and Camille Male for several useful discussions, and the anonymous referee of this 
paper for a number of useful comments and suggestions. The work of T.B. was supported 
by the ANR grant "Granma". I.N. acknowledges financial support from the ANR project 
OSvsQPI 2011 BS01 008 01 and from a CNRS PEPS grant. 



Our starting point is the Marchenko-Pastur theorem [16J. Let us recall that a complex 
Wishart matrix of parameters (N,M) is a random N x N matrix of type W = j^GG*, 
where G is a N x M matrix with independent complex Gaussian Af(0, 1) entries. 

Theorem 1.1. In the limit N, M — > oo, M/N — > t > 0, the law oftW converges to 
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[(Vi-i) 2 ,(Vt+i) 2 ]( x ) d- x 



which is called Marchenko-Pastur law of parameter t > 0. 
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Proof. This follows for instance by checking, by using the Wick formula, that the asymp- 
totic moments of W coincide with those of ir t . We refer here the reader to section 2 below, 
where we present a generalization of this theorem, along with a complete proof. □ 

In order to deal with the block-modified case, we will need a free probability point of 
view on ir t . The idea is that TT t appears naturally via a "free Poisson limit" procedure. 

We recall from [21] that a noncommutative probability space is a pair (A, ip), where A 
is a unital C*-algebra, and if : A — > C is a positive unital trace. The law of a self-adjoint 
element a G A is the probability measure on the spectrum of a (which is a compact subset 
of R) given by J f(x)dfi(x) = (p(f(a)), for any continuous function / : R — y C. 

Two subalgebras B, C C A are called free if <p(. . . biCib i+ iCi + i . . .) = whenever 6, e B 
and qsC satisfy f(bi) = f(ci) = 0. Two elements b,c G A are called free whenever the 
algebras B = (b) and C = (c) that they generate are free. Finally, the free convolution 
operation EH is defined as follows: if fi, v are compactly supported probability measures 
on R, then \x E3 v is the law of b + c, where b, c are free, having laws /i, z/. See [21]. 

With these definitions, we have two conceptual results about n t , as follows. 

Proposition 1.2. We have 7t s+t = K s \BiTt, for any s,t > 0, so that the Marchenko-Pastur 
laws form a semigroup with respect to Voiculescu's free convolution operation. 



Proof. We recall from [21] that the operation EH for real probability measures is linearized 
by the .R-transform, constructed as follows: first, we let f(y) = 1 + m\y + m^y 1 + ... be 
the moment generating function of our measure, so that G(£) = is the Cauchy 

transform; then we set R(y) = K(y) — y^ 1 , where K(y) is such that G(K(y)) = y. 
By Stieltjes inversion, the Cauchy transform of ir t is given by: 

Thus we can compute the i?-transform, by proceeding as follows: 
£G 2 + 1 = (£ + 1 - t)G => Ky 2 + 1 = (K + 1 - t)y 

Ry 2 + y + 1 = (R + 1 - t)y + 1 
R = t/(l-y) 

Now since the expression tj (1 — y) is linear in t, this gives the result. □ 
Theorem 1.3. We have the free Poisson limit formula 

^ ^ \ EBn 



7T t = lim 1 5 + -5i 

n->oo y y n J n 

so that -Kt is called as well 'free Poisson law" of parameter t > 0. 

Proof. This result is once again well-known, see e.g. Speicher [20] or Hiai and Petz [13], 
and is for instance a particular case of Theorem 1.5 below. □ 



BLOCK-MODIFIED WISHART MATRICES 




In what follows we will need some generalizations of the above results. The class of 
"compound free Poisson laws" was introduced by Speicher in [20], studied by Hiai and 
Petz in [13], then further studied in [3], [6]. Since we will be only interested in the "discrete 
case", it is technically convenient to introduce these distributions as follows. 

Definition 1.4. For fi = Y2i=i c $zi with q > and Zi G C, we let 

ix p = law 

where the variables at are free Poisson of parameter Ci, free. 

Observe that we don't necessarily assume fi to be of mass 1. Observe also that we 
don't assume the numbers Zi to be distinct: the fact that we can indeed do so comes from 
Proposition 1.2, which shows that the distribution tt^ is indeed well-defined. 

Theorem 1.5. If fi is real we have the Poisson limit formula 

( ^ \ EBra 
(1 ) 8 + -// I 
V n/ n J 

where c = mass(fi), so that tc p is called "compound free Poisson law" associated to \i. 

Proof. This result is from [5], we present below the idea of the proof. Let p n be the 
measure appearing in the statement, under the convolution sign. We have: 



^(o=(i--H+-x: 

V n) f n 



Now since K Pn {y) = y" x + R Pn {y) = y^ 1 + R/n, where R = R p m n (y), we get: 



c\ 1 1 v— > c 



-y— 



n) y~ l + R/n n ^—^ y~ l + R/n — Zi 
Now multiplying by n/y, rearranging the terms, and letting n — > oo, we get: 

C + yR _ y-\ Ci ^ ^ . . _ CiZi 

1 + yR/ n + U R I n ~yzi " y Zi 

On the other hand, let a be the sum of free Poisson variables in the statement. By 
using the i?-transform formula in the proof of Proposition 1.2, we have: 

s 

R ai (y) = ~, — Rz iai {y) = 1 *' R a (y) x 



1 - y ^-yzi j^i-yzi 

Thus we have indeed the same formula as above, and we are done. □ 
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Finally, we will need the notion of free cumulant, introduced by Speicher in |19j . The 
free cumulants K p (a) of a self-adjoint variable a are the coefficients of its i?-transform: 

oo 

Ra(y) = ^K p+1 (a)y p 

With K n (a) = Ylbtw K #b(a), where the product is over all blocks of it, and # is the size 
of blocks, we have then the following moment-cumulant formula, due to Speicher |19j : 

n£NC(p) 

For concrete applications, it is rather the converse statement that we will use: whenever 
we have a sequence of numbers K p (a) such the above moment formula holds for any k G N, 
it follows that these numbers K p (a) are the free cumulants of a, see [19|. 

The following will be our main tool for detecting compound free Poisson laws. 

Theorem 1.6. If is real then the free cumulants ofir^ are the moments of [i. 

Proof. We write fj, = Yli=i c ^Zi with q > and G R. We know from the proof of 
Theorem 1.5 that the i?-transform of 7r M is: 

s s oo oo / s \ 

w = E = E E «^rv = E E v 

i=l ^ * i=l p=0 p=0 \i=l / 

Now since the free cumulants are the coefficients of R, this gives the result. □ 

2. Wish art matrices 

Consider the embedding NC (p) C S p obtained by "cycling inside each block" . That is, 
each block b = {b±, . . . , bk} with b\ < . . . < bk of a given noncrossing partition n G NC(p) 
produces by definition the cycle (bi . . . bk) of the corresponding permutation n G S p . 

Observe that the number of blocks of 7r G NC (p) corresponds in this way to the number 
of cycles of the corresponding permutation 7r G S p . This number will be denoted \ir\. 

For 7r G NC(p) we denote by ||7r|| the number of blocks of 7r having even size. We will 
need a number of facts on partitions, summarized in the following statement: 

Lemma 2.1. If 7 G NC(p) is the one-block partition, then \tt\ + |7T7 _1 | < p+1, with 
equality iff it G NC(p). In addition, for tt G NC(p) we have: 

(1) |7 _1 vr| = |vt7 _1 | = p + 1 — \tt\ . 

(2) |1| + |7r'7| = |7T7 _1 | + |7r' +1 |, for any I G Z. 

(3) |77r| = |7T7| = 1 1 -TT - 1 1 + 1. 

Proof Observe first that 7 G S p is the full cycle, 7 = (1 . . .p). Let also 1 G NC(p) be the 
p-block partition, so that the associated permutation 1 G S p is the identity. 
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It is known that l(jr) = p — \tt\ is the length of permutations, so that d(7r, or) = li^na^ 1 ) 
is the usual distance on S p . Now the triangle inequality d(l,j) < d{l,Ti) + d(7r, 7) reads 
p — 1 < (p — I vr I ) + (p — |7T7 _1 |), so we have \ir\ + |7T7 _1 | < p + 1, as claimed. For the 
assertion regarding the case where we have equality, see Biane [8]. 

(1) This is clear from the first assertion, or from the well-known fact that 7" 1 vr, 7T7 -1 
have the same cycle structure as the right and left Kreweras complements of tt. 

(2) By using |1| = p and |vr7 _1 | = p + 1 — \tt\, we must prove the following formula: 

1 7T 7| — 1 = |7T I — I vr I 

Observe first that this formula holds for 7r = 7. We will prove this formula by recurrence 
on the number of blocks of tt. Since n is noncrossing, it is enough to check the stability 
of the above formula by "concatenation", tt = (7T1, tt2). Since the right term of the above 
formula is additive with respect to concatenation, we have to prove that the left term is 
additive as well with respect to concatenation, i.e. we have to prove that: 

\(TT 1 ,TT 2 ) l 1 \-l = (\TT l lll \-l) + (\TT l 2l2 \-l) 

With o"! = tt\ and a 2 = tt 1 2 , the formula to be proved becomes: 

1(0-1,0-2)7! - 1 = (kiTil - 1) + (I0272I - 1) 

In order to prove this latter formula, we use the following key identity: 

(71,72) = j(k,p) 

More precisely, by using this identity, we have as claimed: 

ki7i| + |o 2 7 2 | = 1(01,0-2) (7i>72)| = l(oi,o 2 )7(A;,p)| = I (o-i, 0-2)7! + 1 

Here at right we have used the general fact that |vr(a, b)\ — \tt\ + 1 when a, b are in the 
same cycle of tt, and |vr(a, b)\ = |vr| — 1 if a, b are in different cycles of tt, and the fact that, 
in the above situation, k,p must be in the same cycle of (01, 0-2)7. 

(3) This is the combinatorial lemma from our previous paper [5] , which follows as well 
from (2). Indeed, at I — 1 the formula in (2) becomes |1| + 1 7r^y | = | Tr^y" 1 1 + |7i~ 2 |, so: 

I vr'y I = 1 7r^y ~ 1 1 + |7T 2 | — |1| 

= (p + 1- \7t\) + (\<k\ + \\7t\\)-P 

= IKII + i 

Together with the fact that |7T7| = |7vr|, this gives the result. □ 

Consider now a complex Wishart matrix of parameters (dn, dm), W = -^GG*, where 
G is a dn x dm matrix with independent complex A^(0, 1) entries. As explained in the 
introduction, we are interested in the study of the "block-modified" versions of W , which 
are defined as follows: 

Definition 2.2. Associated to any Wishart matrix W of parameters (dn,dm) and any 
linear map <p : M n (C) —> M n (C) is the "block-modified" matrix W = (id <8> <p)W . 
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• lfi(A) 3 




^« 3 — • A ^ 

A r I 1 



W 




(b) 



Figure 1. (a) The action of the map (p through A, and (b) the block- 
modified Wishart matrix W as the action of A on W. 

In what follows, the idea will be to relate the limiting d — > oo eigenvalue distribution 
of W to the usual eigenvalue distribution of <p, viewed as a square matrix. 
In order to view (p as a square matrix, we use the following correspondence: 



Definition 2.3. Associated to any linear map (p : M n (C) — > M n (C) is t/ie square matrix 
A G M n (C) <g> M n (C) piuen 6y 

= (Tr <g> id)[(t <g> id) A • (A <g> 1)] 
where Tr and t are the usual trace and transposition of the n x n matrices. 

Observe that the correspondence <p — > A is bijective. We will use as well the correspon- 
dence in the other sense: every time we will have a matrix A e M n (C) <S> M n (C), we could 
speak about the associated linear map ip : M n (C) — > M n (C), via the above formula. This 
correspondence, when restricted to subclasses of (completely) positive maps, is known in 
quantum information theory as the Choi-Jamiolkowski isomorphism. See [7]. 

In order to get now more insight into the multiplicativity condition, we use the tensor 
planar algebra of Jones [H], or, equivalently, the diagrammatic formalism in [10J. We 
represent the matrices A e M n (C) <8> M n (C) as vertical boxes with 2 left legs and 2 right 
legs, so that the correspondence A — > <p is the one described by Figure 1 (a) , where we 
write A r = (t <8> id)A for the partially transposed map. 

It will be often convenient to write A in usual matrix notation: 

A = ^2 ^ab,cde ac <S> e bd 

abed 

With this convention, the associated linear map <p is given by: 

<p{A)bd = ^ A abtCd A ac 



Observe also that, in terms of the above square matrix A, we have (see also Figure 1(b) 
for the diagram): 

W = (id ® Tr ® id){(l ® (t®id)A)(W® 1)] 
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Finally, observe that, in matrix notation, the entries of W are given by: 

Wi a j b = A eaJb W ie jf 

Given a complex linear functional E : A — > C and a permutation f 6 S 1 ,, we set 
E,r(a) = XlbeTr E(a* 6 ) for any a G A, where the sum is over all the cycles of n. We will use 
this notation throughout the reminder of the paper, often with E being the usual trace of 
matrices, Tr : Mjy(C) — > C, or the normalized trace of matrices, tr : Mjv(C) — > C. 

Theorem 2.4. We have the asymptotic moment formula: 

\im(Eotr)((mW) p ) = V (mn) |7r| tr (7r>7) (A) 

ci— 5>oo ' J 

tt€NC( P ) 

Proof. According to the above formula for the entries of W, we have: 
tr(Wn = (dn)- 1 II 

= (dn)- 1 II A — ; -'- ir '~'-- '^ 

%f CLf Gf" f'p & 

= (dnY l (dm)- p ^ JJ K s a s j s a s+1 G isea>jsbs G is+lfa>jshs 
The average of the general term can be computed by the Wick rule: 

E ^11^''-' •; '' ( '''- - : / .;-'• ^ = #i n E Sp\in(s) = is+1, eir(s) = fs, jn(s) = js, K(s) = b s } 

Let us look now at the above sum. The i, j, b indices range over sets having respectively 
d, d, m elements, and they have to be constant under the action of 7T7 -1 , tt, it. Thus when 
summing over these i, j, b indices we simply obtain a d' 71 " 7 'rf'^'m' 71 "' factor, so we get: 

7rSS p a r e r s 

= rr l rrT v d W+K7 ~ lhp_1 m w Tr (7rj7) (A) 
= mT v d kl+k7 ~ 1| ~ p ~ 1 (mn) |7r| tr ( ^ i7) (A) 

By Lemma 2.1, with d — > oo the sum restricts over ir e NC(p), and we are done. □ 

One can also prove the result using the graphical Wick technique developed in [IT] . 
There is an obvious relation between the above result and the combinatorics of the com- 
pound free Poisson laws. More precisely, we have the following similar result: 
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Theorem 2.5. If A = A* then the moments ofn mnp , with p = law(A), are: 

M p = (rnnptr M (A) 

n£NC(p) 

Proof. We know from Theorem 1.6 that the free cumulants of 7T mnp are the moments of 
mnp, given by k p = mn ■ tr(A p ) = mn ■ ir( 7j7 )(A). Together with Speicher's moment- 
cumulant formula, explained in section 1 above, this gives the result. □ 

In what follows we will try to exploit the obvious similarity between Theorem 2.4 and 
Theorem 2.5. We will split our study into 2 parts: the unital case will be investigated in 
section 3, and the general, non-unital case will be investigated in section 4. 

3. The unital case 

A linear map <p : M n (C) — > M n (C) is called self-adjoint if A = A* implies <f(A) = <p(A)*. 
This is the same as asking for the corresponding matrix A to be self-adjoint. 

We use the notation E p from the previous section, in the case where E = tr is the trace, 
and p = (tt, a) is a product of permutations, via the embedding S p x S q C S pq . 

Theorem 3.1. Assume that ip is self- adjoint. Iftr^^A) = tr( 7rj7r )(A) for any p e N and 
any n G NC{p), then mW ~ TT mnp , with p = law (A). 

Proof. We know from Theorem 2.4 and Theorem 2.5 that the formula mW ~ ^ m np is 
equivalent to the following formula, which should be valid for any p G N: 

(mn)Hfr M (A)= ^ (mnptr^A) 

Now since in the case tr^ n ^(A) = tr^ n ^(A) this formula holds, we are done. □ 

The point now is that, in all the examples that we have, the above result applies only 
when if is unital modulo scalars. The non-unital case, which requires a more subtle 
combination of Theorem 2.4 and Theorem 2.5, will be discussed in the next section. 

For a matrix A G M n (C) we denote by A 5 G M n (C) its diagonal. 

Theorem 3.2. The formula tr^^(A) = tr^^A) holds for any map which is unital 
modulo scalars, in the sense that ip(l) = cl with c G C — {0}, of the following form: 

(1) ip(A) = Tr{BA)C. 

(2) tp(A) = BAC. 

(3) tp(A) = BA l C. 

(4) tp(A) = B(EAF) S C. 

In addition, the set of above maps is stable by composition and inversion, and the formula 
tr^^) (A) = tr^^) (A) is stable by taking tensor products, and multiplying by scalars. 



BLOCK-MODIFIED WISH ART MATRICES 11 

Proof. If we write A = J2 a bcd Kb,cd e ac® e bd , then (p(A) bd = J2ac A ab,cdA ac . This shows that 
in the cases (1-3) we have A ab)Cd = B ca C bd , B ba C cd , B bc C ad respectively. 
Let us compute the numbers Tr (7r , CT) (A) = E„ r e r EL A e s a s , Ms)Ms) : 

Q-T' Cf S 

Tr M (A 2 ) = EII fl ^.)M.) = Vi(BC) 
Tr (7r , CT) (A 3 ) = ^n^w^^Tr^^C) 

In the case (4) now, we have A abjCd = Y, x B bx E xa F cx C xd , and we get: 

^ r (vr,o-)(A ) = ^ ] J^J B asXs E Xses F e ^^ Xa C Xsa<T ^ = ^ ] | | (^ < -^)x s a tT ( S ) (EF) Xn(s)Xs 

a r e r x r s x r s 

Assume now that the maps ip are unital modulo scalars, with C = cl in (1), with 
BC = cl in (2,3), and with {CB){EF) S = cl in (4). We get that: 

tr^A 1 ) = tr w (B)tr a (C) = cHr^B) 

tr M (A 2 ) = nl^'l-l-l-H^^^c) =n l-- 1 l-kl-kl c P 

tr ( ^ i(T) (A 3 ) = n^-^-^kr^BC) = n M-M-M c P 

tr M (A 4 ) = £ (CB) xt (EF) xx = n^-^<f 

x<kcr(nA<r) 

We claim that we have tr^ tJ ^(A) = tr^^{A) in all cases. Indeed, in case (1) this is 
clear, and in (2-4) this follows from the following equalities for the exponents of n: 

(j) + 1 — 1 7T | ) — 1 7T | — 1 = p — 1 7T | — 1 7T I 

(||7r|| + 1) - |vr| - 1 = + \\ir\\) - \n\ - \n\ 

1 — 1 7T | — 1 = | -7T | — | VI" | — |/T| 

Regarding now the stability properties, the inverses of the maps (1-4) are also of type 
(1-4). Also, the maps (1-4) are stable by composition, the formula for i o j being: 

i\j 12 3 4 
11111 

2 12 3 4 

3 13 2 4 

4 14 4 4 
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For the tensor product assertion now, assume that ip = f 1 <8> <fi 2 - We have then A 
A 1 <g> A 2 , so in matrix notation A aab p :Cyd s = A abcd A 2 a p 1&1 and we get: 



y ttal e a At 



= Tr ( ^ (T) (A 1 )Tr (7r , CT) (A 2 ) 

It follows that tr( 7r o .)(A) = tr( 7r)0 -)(A 1 )tr( 7r o -)(A 2 ), and we are done. Finally, the assertion 
regarding the multiplication by scalars is clear. □ 

As a first consequence, we obtain some previously known results, from [TT], [IB], [5]: 

Proposition 3.3. We have the following results: 

(1) t(id g) tr(.)l)H / ~ ir t , where t = mn. 

(2) tW ~ 7T t; where t = m/n. 

(3) m(id <S> t)W ~ law(a — (5), where a, /3 are free Poisson (m(n ± l)/2) ; /ree. 

(4) m(zrf® ~ vr m . 

Proof. Our claim is that (1-4) above correspond via Theorem 3.1 to the assertions (1-4) 
in Theorem 3.2, at B — C — 1. Indeed, let us see what happens in this case: 

(1) Here we have <p(A) = Tr(A)l, hence W = (id ® Tr(.)l)iy. Also, we have A abjCd = 
S ca 5 bd , hence A = ^ a6 e ao <g)e 66 is the identity matrix: A = 1. Thus p = 5 1 , so n mnp = n mn , 
so Theorem 3.2 says at B = C = 1 that we have mW ~ vr mn , as claimed. 

(2) Here we have <p(A) = A, so W = W. Also, we have A ab ^ cd = 5 ab 5 cd , hence A = 
Ylac e ac ® e aci so we have A = nP, where P is the rank one projection on the vector 
£ a e a ® e a e C n ® C\ Thus p = ^6 + ^5 n , so mnp = ^=^5 + *5 n , so vr mnp = 
law(na), where a is free Poisson (m/n). Thus Theorem 3.2 says at B — C — 1 that we 
have miy ~ law(na), and by dividing by n we get (m/n)W ~ law(a), as claimed. 

(3) Here we have <p(A) = A f , so If = (id <S> t)W . Also, we have A abfid = 5 bc 5 ad , 
so A = J2ac e ac ® e ca is the flip: A(e c (8> e a ) = e a <8> e c . Thus p = ^6-1 + so 
mnp = m ^ 1 - > (5_ 1 + m (^ +1 ) ^ 1; go 7r mnp = law(a — /3), where a, /3 are as in the statement. 
Thus Theorem 3.2 says at B — C — 1 that we have mlf ~ /aw (a — /3), as claimed. 

(4) Here we have <p(A) = A 5 , so W = (id® (.) 5 )W . Also, we have A abcd = b~ a ,b,c,di hence 
A = e aa <8> e aa is the orthogonal projection on span(e a ® e a ) C C n <E> C n . Thus we have 
p = ^E^So + so mnp = m(n — 1)5 + m8\, so 7r mnp = n m . Thus Theorem 3.2 says at 
B = C = 1 that we have rnW ~ ir m , as claimed. □ 

We have as well the following generalization, obtained by taking a tensor product: 

Proposition 3.4. For a linear map of type ip = id g) (tr(.)l ®> id ® t ®> (.) 5 ) we have 
tW ~ law(a — j3), where a, (3 are free Poisson (t(n 3 ± l)/2), free, with t = mn\/n 2 . 
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Proof. We have A = (n2/rai)(l £g> P <g> E £§> Q), where P, £, Q are the matrices appearing 
in the above proof, namely: P is a rank 1 projection in n\ dimensions, £ is the flip in n\ 
dimensions, and Q is a rank projection in n\ dimensions. Since P ® Q is a rank 
projection in (^n^) 2 dimensions, the matrix P (g) X (g) Q follows the following law: 

_ n|n4 - 1 ^ n 3 -l ^ n 3 + 1 ^ 

Now since tensoring with id doesn't change the law, it follows that p is the measure 
obtained from rj by replacing the ±1 atoms by ±n 2 /ni atoms, so we have: 

mn\n 3 {n\n^ — 1) mni(n 3 — 1) mni(n 3 + 1) 
mnp = S + S_ m/m + 5 m/ni 

By using Theorem 3.1 and Theorem 3.2 it follows that mW ~ law((ri2/ni)(a — (3)), 
where a, (3 are as in the statement, and by multiplying by n 1 /n 2 we obtain the result. □ 

4. The main result 

In this section we investigate the general, non-unital case. Let us first work out a 
generalization of Theorem 3.1, which will prove to be well-adapted to this case. 

Theorem 4.1. Assume that ip is self-adjoint. If for any p G N and tt G NC(p) we have 

*T(i ) i)(A)tr( w , 7 )(A) = £r ( i )7r7 -i)(A)tr (7ri7r )(A) 

then 5mW ~ vr mnp IE v , with p = law(A), v = law(D), 5 = tr(D), where D = <p(l). 

Proof. The entries of D = </?(l) are given by D bd = ^ a A a6 arf , so we get: 

tr a {D) = J2H D ^ (r) = A ^r,ae a(r) = *r (v) (A) 



In particular at o = 1 we obtain 5 P = trhi)(A), so Theorem 2.4 gives: 
lim (E o tr)((5mW) p ) = V (mn) w tr (lil) (A)tr K7) (A) 



d— >0O 

TreNC(p) 



We use the following formula from [T7j for of a multiplicative free convolution: 

tt£NC(p) 6Stt 

We know from Theorem 1.6 that the free cumulants of ir mnp are the moments of mnp, 



7T 



mnp 



given by k p = mn ■ tr(A p ) = ran ■ £r( 7>7 )(A). Thus the abstract free cumulants of 
are the numbers = (mn)' 7r 'tr( 7rj7r )(A). Together with the formula tr a (D) = £r(i i<7 )(A) 
above, the above general formula from [ 1 Tj shows that the moments of TT mnp E v are: 

M p = J2 (mn) |7r| tr ( i i7r7 -i)(A)tr (7r)7r) (A) 

n£NC(p) 
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We conclude that the formula SmW ~ ir mnp M v is equivalent to: 

(mn)Mtr (lil) (A)tr (7ri7) (A) = (rnn^kr^-i^tr^A) 

tt£NC(p) n£NC(p) 

With this equivalence in hand, the assertion in the statement is now clear. □ 

The above result suggests the following technical definition: 
Definition 4.2. A map (p and the corresponding matrix A are called "multiplicative" if 

tr(i ) i)(A)ir( w>7 )(A) = tr ( i i7r7 -i)(A)tr (7ri7r) (A) 
for any pGN and ir G NC (p) . 

Observe that in this definition one can replace the normalized traces tr by the unnor- 
malized traces Tr: this follows indeed from Tr( 7rjCr )(A) = n^^Hr^^A). 

As already mentioned, this definition is a technical one, coming straight from Theorem 
4.1, which in turn comes from abstract algebraic manipulations. There is no simpler 
formulation of it, but in what follows we will try to have some understanding of it. 

Let us extend now Theorem 3.2 above. We have the following result here: 

Theorem 4.3. The following types of maps and matrices are multiplicative: 

(1) (p(A) = Tr{BA)C, or A = B®C, in the case C = cl. 

(2) ip(A) = BAC, or A = \B)(C\, for any B,C. 

(3) tp(A) = BA l C, or A = SWAP BC; in the case BC = cl. 

(4) (p(A) = xA 5 , or A = Center^,, in the case x = cl. 

In addition, the set of multiplicative maps and matrices is stable by tensor products. 

Proof. (1,3,4) follow from Theorem 3.2, and the last assertion follows from the proof of 
Theorem 3.2. For proving (2) we use a formula found in the proof of Theorem 3.2: 

Tr M (K) = Tr na -,(BC) 

In terms of normalized traces, and with D = BC, we get: 

tr M (A) = n^- 1 \-^tr 7Ta - 1 (D) 

In particular we have the following formulae: 

tn,i(A) = n p - p - p t ri (D) =n- p tn{D) 
tr^(A) = #i-kl)-kl-i tvi p) = nV-^kr^-iiD) 
tr^ 7 -i(A) = n^ +1 -^~ p ~^ +1 -^hr^-i(D) = n- p tr^-i(D) 
tr ntW (A) = rf-W-Wtr^D) = rp-Wtr^D) 
By multiplying we obtain the formula in the statement. □ 
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Figure 2. Special types of matrices A whose mixed moments factorize properly. 

Using the graphical notation from [10], the diagrams for the maps of type (1-4) in 
Theorem 4.3 are some very simple ones, namely those in Figure [2j 

So, let us try to combine the diagrams in Figure 2 in a "planar" way. We call "gener- 
alized Center c diagram" the diagram A having k left legs and I right legs, given by: 

The following planar-categorical definition is inspired from [H], [TO] : 

Definition 4.4. The "strings and beads" operadV is defined as follows: 

(1) The elements are the matrices A G (M rtl (C) (g) . . . (g) M nk {G))® 2 , represented by 
vertical boxes with 2k left legs and 2k right legs. 

(2) The legs are colored by the values of the corresponding numbers, and the tensor 
product and composition operations have to match colors. 

(3) Inside the box we have "strings" joining the 2k + 2k legs, i.e. we have a partition 
of 2k + 2k elements into even blocks tt G P even (2k,2k) , with colors matching. 

(4) Each string can be decorated with "beads", i.e. with usual matrices A G M n (C), 
where n is the color of the string. 

(5) The "multi-leg" strings, representing blocks of tt G P even (2k,2k) having size > 4, 
have at the multi- crossing a "small bead", of generalized Center c type. 

In other words, the matrices A G (M ni (C) <g) . . . <g) M nfe (C))® 2 will be represented as 
diagrams having the left and right sequences of 2k points marked n\, . . . , n^, n^, . . . , n±, 
from top to bottom. In order to view such a matrix as A G M ni ... nfe (C) ® M rii ,... ril (C), we 
will simply "compact the blocks" , i.e. we will view A as a diagram between 2 left points 
and 2 right points, with both 2-series labeled ni . . . n^, Uk ■ ■ ■ n\, from top to bottom. 
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We know from Theorem 4.3 that the 4 types of diagrams appearing there, as well as 
their tensor products, are multiplicative in the sense of Definition 4.2. The challenging 
question is to determine all the elements of V which are multiplicative. 

Let us call "through" strings of A G V the strings of the corresponding partition, joining 
left and right points. With this convention, we have the following result: 

Proposition 4.5. Assume that A <EV is self-adjoint, and that its through strings are not 
decorated with beads. Then A is multiplicative. 

Proof. The element A consists by definition of a partition x G P even (2k,2k), decorated 
with beads. Since A is self-adjoint, x is "symmetric", in the sense that it is invariant 
under the reflection with respect to the vertical middle axis. We decompose x = Uj e j&j, 
where hi are blocks or pairs of blocks of x, chosen symmetric, and with |/| maximal. 

For any such "symmetric block" fej we denote by A* G V the element obtained from A 
by keeping 6j with its beads, and by completing with (undecorated) horizontal strings. 

It follows from definitions that we have: 



This formula shows that the multiplicativity condition splits over the symmetric blocks. 
Now for the symmetric blocks consisting of pairings, the result follows from Theorem 4.3 
above. So, it remains to prove that the symmetrization A of any generalized Center c type 
partition A is multiplicative. By adding an extra string if needed we can assume that A 
has an even number of legs at left and at right, and we have two cases here: 

(1) A = A. Here the formula of A is as follows, with m G {2k — 1, 2k}: 



^)(A)=n Tr (^)(^) 



■iik,n—3ik 



= c • 5i 



l2k,Jm+l---32k 



With these notations, 



the mixed trace is given by: 





c p . 7l l 7rA °1 . n ( 2fe - m )M 



Now^ since both nl 7rAcr l and are multiplicative, this gives the result. 
(2) A^A. Here the formula of A is as follows, with m G {2k — 1, 2k}: 

1V h...i2k,h—32k ~ C ' h,--;is,ik+l,--;ik+t,js + l, — ,jk:jk+t+l — ,jm 

r 

°i s+ i,...,i k ,i k+t+1 ,...,i m ,ji,...,j s ,j k+1 ...,j k+t 
&im+l---i2k,im+l—32k 
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With these notations, the mixed trace is given by: 

TV(7r,<r)(A) = /] TT C <^ il ii i< Ml) Ml) MD 

<^-i ji ji Mi) Mi) Mi) Mi) 

By looking at the indices we obtain Tr( 7ri(7 )(A) = c 2 • n 2 ' 71 ^ "' ■ 7i( 2fe_m )l°1 ) and since both 
the quantities n' 7rA<T ' and n'°"' are multiplicative, we are done. □ 

We can state and prove now our main result: 

Theorem 4.6. Let W be a complex Wishart matrix of parameters (dn,dm), and let 
ip : M n (C) M n (C) 6e a linear map coming from "strings and beads" diagram A e V, 
which is symmetric and has undecorated through strings. Then the d — > oo asymptotic 
eigenvalue distribution of the block-modified Wishart matrix W = (id ® ip)W is given by 
5mW ~ iimnp v, with p = law(A), v = law(D), 5 = tr(D), where D = <p(l). 

Proof. We know from Proposition 4.5 that A is multiplicative in the sense of Definition 
4.2. Thus Theorem 4.1 applies, and gives the result. □ 

Observe that Theorem 4.6 doesn't fully cover Theorem 4.3. The problem is that, when 
trying to fully combine the diagrams in Theorem 4.3, our formula SmW ~ 7r mnp IEz/ seems 
to need a kind of substantial upgrade, with more quantities involved on the right. This 
kind of upgrade is also suggested by the various questions raised in the introduction, but 
finding it is of course a quite difficult problem, that we would like to raise here. 
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